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ABSTRACT 

Aims. Recent studies showed that time-distance inversions for flows start to be dominated by a random noise at a depth of only a few 
Mm. It was proposed that the ensemble averaging might be a solution to learn about the structure of the convective flows, e.g., about 
the depth structure of supergranulation. 

Methods. Time-distance inversion is applied to the statistical sample of ~ 10 4 supergranules, which allows to regularise weakly 
about the random-noise term of the inversion cost function and hence to have a much better localisation in space. We compare these 
inversions at four depths (1.9, 2.9, 4.3, and 6.2 Mm) when using different spatio-temporal filtering schemes in order to gain confidence 
about these inferences. 

Results. The flows inferred by using different spatio-temporal filtering schemes are different (even by the sign) even-though the 
formal averaging kernels and the random-noise levels are very similar. The inverted flows alterates its sign several times with depth. 
It is suggested that this is due to the inaccuracies in the forward problem that are possibly amplified by the inversion. It is possible 
that also other time-distance inversions are affected by this issue. 
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1. Motivation 

Helioseismology is considered to be a state-of-the-art method 
of solar research to learn about the solar interior (see review 
by |Gizon et al.||2010| }. Among various helioseismic methods, 
the local helioseismology plays an important role in studying 
the spatially localised details. The time-distance helioseismol¬ 
ogy ( [Duvall et al.|1993| is one of the approaches. It consists of 
tools for measuring and interpretation of travel times of solar 
waves. Since its formulation, it was used extensively to infer es¬ 
pecially information about flows and sound-speed anomalies in 
upper layers of solar convection zone. 

The method assumes that the travel time r a is bound to the 
plasma anomalies Sq a with respect to the background by a linear 
equation 


5r a (r) = 



Kp(y - r)6qp(r f ,z) + n a (r). 


( 1 ) 


where K% is the sensitivity kernel coming from the forward mod¬ 
elling, r = (x,y) is a horizontal position vector in a Cartesian 
coordinate system (the remaining vertical component will be de¬ 
noted as z) and n a is the travel-time random noise. Index a de¬ 
notes the selection of the waves (the combination of the spatio- 
temporal filter, spatial averaging, and the distance between the 
travel-time measurement points) and greek indices a or (3 select 
the perturbation 6q (flow components, sound speed perturbation, 
...). 

One of the main goals of time-distance helioseismology is to 
infer the structure of convective flows in the near-surface layers 


of solar convection zone from surface measurements of wave 
travel times. Formally it means that we replace dqp in 0 by vp, 
where fi = (x, y, z), and search for it. This can only be achieved 
by inverse modelling. 

To solve the inverse problem, one has to set up the cost func¬ 
tion with various terms, which can be cast to a system of linear 
equations. How to perform such task was in details described and 
discussed elsewhere. The resulting flow estimate v inv is given by 


v7fa>;zb) 



( Kp(r - r 0 , z; Zo)v,s(x) d 2 r dz 


+ E w °( ri ~ r ° ; Zo )n a (ri), 

i,a 


( 2 ) 


where w a a are the inversion weights to be determined and *Kp in¬ 
dicates the component of the inversion averaging kernel (a linear 
combination of the sensitivity kernels with inversion weights). 
The averaging kernel describes the smoothing of the real solar 
convective flows and is constructed to peak around a target depth 
Zo- The second term on the right-hand side of ^ indicates the re¬ 
alisation of the random noise. 

For the inversion method which is t he bas e of this study, the 


reader is referred to Jackiewicz et al. (2012) and 


Svanda et al. 


( |2011| ). The cost function contains four terms balanced by three 
trade-off parameters: the quality of the fit (the misfit) of the av¬ 
eraging kernel to a user-defined target function localised in the 
Sun, the level of the random noise, the level of the pollution of 
the inverted flow component by other components (the cross¬ 
talk ;), and the term which ensures that the resulting inversion 
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weights are spatially localised (which is a requirement in order 
to fulfil the mathematical assumptions). 

One of the biggest issues in time-distance helioseismology is 
the presence of the realisation noise. Pressure and surface grav¬ 
ity waves are randomly excited by the vigorous convection and 
hence their statistics inherently contains a large random com¬ 
ponent. This acts as a random noise in any of the helioseismic 
observables and travel times of the waves are no exception. Var¬ 
ious helioseismic methods deal with the random noise differ¬ 
ently. Sometimes the noise is ignored (e.g. |Kosovichev & Duvall| 
11997k sometimes the variance of the travel times is considered 
(e.g. Zhao et al.|2001| ) and in the ideal cas e, the full travel-time 
noise covariance matrix is considered (e.g.|Couvidat et al.|2005 ' 


Jackiewicz et al. 2008 Svanda et al.|2011 ). 

The precise knowledge of the noise covariance matrix is es¬ 
sential, as this error propagates through the inversion procedure 
and translates into the realisation of the random error in the 
inverted maps. The determination of the full noise-covariance 
matrix is not straightforward. There are two approaches being 
used. Either the covariance matrix is read out directly from the 
data by measuring a large set of travel-time maps, or it is de¬ 
rived from the model using a Monte-Carlo-like approach. Both 
approaches have their own advantages and disadvantages. The 
data-driven approach would be probably ideal, however there 
might be time-varying systematic errors in the travel-time mea¬ 
surements, which will affect also the estimated error level of the 
inverted quantity. The model-driven methods relies on how pre¬ 
cisely does the used model agree with the observed power spec¬ 
trum of the waves. So it would be ideal if one could minimise 
the effect imprecise knowledge of travel-time noise to the knowl¬ 
edge of solar flows. 

Moreover, several studies ( Woodard 2007; Jackiewicz et al. 
2008) showed that the inversions for snapshot of the solar flows 
averaged over a few hours starts to be dominated by the random- 
noise com ponent already at very shallow depths. 


In 


Svanda et al. 


(2011) we suggested the use of the ensem¬ 


ble averaging approach. The ensemble averaging uses a strong 
assumption that the random-noise realisations in the individual 
representatives of the ensemble are independent, and hence in 
the average the random-noise level scales as 1 / yfN, where N is 
the number of representatives. When N is large enough, it allows 
to relax the noise term in the inversion cost function and thus to 
decrease the effect of the possible inaccurate knowledge of the 
noise covariance matrix. The ensemble averaging approach cer¬ 
tainly is not useful for studying the snapshots of the flows, how¬ 
ever when investigating the set of the representatives of the same 
phenomenon (supergranules, sunspots,...), it seems extremely 
powerful. In recent years it seems to be a standard method of 
helioseismic research (e.g. |Duvall et al.||20061 |Duvall & Birch| 


2010 §vanda|2012 Birch et al.|2013 Svanda et ai"|2014). 


2. Inversion 

The results presented in this note were obtained using the real 
data coming from the HMI archive. Using the standard tracking 
and mapping pipeline (created and maintained within the Ger¬ 
man Science Center for SDO by H. Schunker and R. Burston), 
24-hours long Dopplergram datacubes were tracked on a daily 
basis, from 8 May 2010 to 12 July 2010. These datacubes cov¬ 
ered the central part of solar disc roughly having 60 degrees on 
a side in the Postel projection with a cadence of 45 s. Each dat- 
acube was processed in a standard way in a travel-time mea¬ 
surement pipeline. From each frame of the datacube, the mean 
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image capturing mostly the pattern of supergranulation was re¬ 
moved. The datacubes then underwent the spatio-temporal filter¬ 
ing for a set of filters (/ to p 4 ridge filters and eleven standard 
phase-speed filters) to retain only the waves with desired proper¬ 
ties. The travel times were then measured from the filtered dat¬ 
acubes using the linearised |Gizon & Birch| ( |2004| ) approach for 
a set of distances for each of the spatio-temporal filters using 
the centre-to-annulus and centre-to-quadrant averaging schemes. 
Travel-time maps of the waves sensitive towards the surface (the 
/ mode and first four phase-speed filters) in a centre-to-annulus 
geometry show clearly the pattern of divergence centres corre¬ 
sponding to supergranulation. 

These travel-time maps were inverted in an inversion 
pipeline. The inversion code is written in Matlab language and 
described in details in Svanda et al. IIzoIB- it utilises the Born- 
approximation sensitivity kernels (Birch & Gizon 2007^ and 
full travel-time covariance matrix measured directly from a large 
set of travel-time maps. The inversion implements the Multi¬ 
channel Subtractive-OLA approach ( Jackiewicz et al.|2012| ) with 
additional terms of the cross-talk minimisation and inversion 
weights localisation. The cost function is regularised strongly 
about these terms, as they both are possible sources of bi¬ 
ases. Both the travel-time and inversion pipelines were validated 
against the surface measurements ( Svanda et al.]2013). 

Three different combinations of sensitivity kernels were 
used. The first set combined all ridge-filtered kernels (/ to p 4 
with a set of annulus radii in the range of 5 to 20 pixels with a 
step of 1 px, hence 240 independent kernels), second all phase- 
speed-filtered kernels (eleven standard phase-speed filters with 
five distances for each of them as tabulated in Table 1 of lCouvi-1 
|dat et al.| ( |2006[ ), together 165 kernels) and finally the inversion 
that combined all these kernels at once (405 kernels together). 


recently (Svanda! 2013, 

DeGrave et al. 

2014b a) and the tests 

by Svanda(2013; 

) showed that for a testing depth they provided 


flow estimates that were highly comparable. These tests however 
were performed in a different regime, when the cost function was 
regularised strongly about the random-noise term. 

A different regime was used in this study by running the 
inversion suitable for the ensemble averaging approach. In this 
case, the requirement on the random-noise term of the inversion 
cost function may be relaxed, which allows the inversion to find 
the solution with a smaller misfit. In another words, the inversion 
results in the averaging kernel that is better localised in space and 
has fewer sidelobes usually making the interpretation difficult. 

The testing sample consisted of 17 474 individual supergran- 
ular cells, identified by a watershed algorithm in the centre-to- 
annulus travel-time maps. The segmentation algorithm and its 
implementation is described elsewhere (Svanda et al. 2014). The 
resulting flow maps of all three flow components were aligned 
about the centres of all detected supergranules and averaged, 
however for the following discussion only the v x component 
(the component in the direction of solar rotation) was used. The 


1 All point-to-point sensitivity kernels were computed by a code kc3 
written by Aaron Birch, which uses the normal mode summation ap¬ 
proach with considered modes upto radial order of 8, quadrupole source 
at the depth of 100 km, and the observation height of 300 km. The cor¬ 
relation time was chosen to be 48 s. A very fine grid in both the wave- 
number (6 x 10 -5 Mm -1 ) and frequency (43 pHz) space was chosen 
for a initial computation, however the results do not seem to be ex¬ 
tremely sensitive to the grid selection. The point-to-annulus and point- 
to-quadrant kernels were computed subsequently utilising spatial aver¬ 
aging. 
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Fig. 1. The inversions for the horizontal flows in an average supergranule at depths 1.9 and 2.9 Mm, together with the display of the corresponding 
averaging kernels. The left column represents the inversion utilising combined ridge+phase-speed filtering, the middle column is for the ridge 
filters only and the right column for phase-speed filters only. 


claims are exactly the same for the remaining horizontal com¬ 
ponent v y and the issues to be revealed get even worse for a 
weak vertical v z component, for which the random-noise con¬ 
straint was relaxed beyond the reasonable signal-to-noise ratio. 

3. Interpretation issues 

The aim of this study was to test the credibility of the time- 
distance inversions for convective flows. In the misfit-dominated 
regime the main burden of the inversion quality is carried by the 


sensitivity kernels which come from forward modelling. Only a 
slight regularisation is performed by the random-noise term. A 
strong regularisation about the cross-talk and weights localisa¬ 
tion terms should guarantee the minimisation of biases. 

All three different travel-time filtering schemes in the inver¬ 
sion provide similar averaging kernels and similar random-noise 
levels. Equation © implies that in that case also the estimate of 
the inverted flow should be similar. Such claim was verified by 
varying the trade-off parameters for the inversions utilising the 
same set of measurements. 
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Fig. 2. Same as Fig. [I] but for the depths of 4.3 and 6.2 Mm. 
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The inverted flow maps alongside with the corresponding av¬ 
eraging kernels and noise levels are presented in Figs.[l]and|2] It 
seems that the above mentioned implication does not hold when 
different sets of measurements are used in inversions. For each 
investigated depth (1.9, 2.9, 4.3, and 6.2 Mm) one can obtain 
different answers (even by a sign) just by selection of the dif¬ 
ferent set of measurements, when the random-noise levels are 
comparable (within the factor of two, given the magnitudes of 
the inverted flows with very large signal-to-noise ratios) and the 
averaging kernels are very similar. It has to be noted that due to 
the strong regularisation of the cost function about the cross-talk 
term, the components of the averaging kernels that are not in the 
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direction of the inversion (equation ([2]), first term on the right- 
hand side for a ± J3) are negligible (not displayed). The second 
apparent issue is that within the inversions based on the same set 
of measurement, the sign of the flow reverses repeatedly as one 
goes deeper. This does not seem to be supported from the side of 
the theory of solar convection. 

The inversions for statistical ensembles are sort of extreme. 
In this extreme regime we however see clearly that they are not 
robust. There is no reasons not to believe that in the case when 
the noise-regularisation term is stronger and the resulting flow 
maps are seemingly realistic, the problems persist to an unknown 
(possibly smaller) extent. For instance the number of reversal 


































Michal Svanda : Issues with time-distance inversions for supergranular flows (RN) 


Inversion 

Ptr-o 

RMS(vf ern ) 
[m s -1 ] 

expected noise 
[ms -1 ] 

derived noise 
[ms -1 ] 

1.9 Mm, combined 

0.99 

70 

14 

11 

2.9 Mm, combined 

0.99 

59 

9 

8 

4.3 Mm, combined 

0.99 

49 

7 

6 

6.2 Mm, combined 

0.99 

45 

6 

5 

2.9 Mm, ridge 

0.99 

71 

14 

12 

2.9 Mm, ridge 

0.99 

59 

12 

9 

4.3 Mm, ridge 

0.99 

49 

10 

8 

6.2 Mm, ridge 

0.99 

45 

8 

6 

1.9 Mm, phase-speed 

0.98 

70 

16 

14 

2.9 Mm, phase-speed 

0.98 

59 

15 

13 

4.3 Mm, phase-speed 

0.97 

49 

15 

12 

6.2 Mm, phase-speed 

0.97 

44 

15 

11 


Table 1 . Tests of the inversion using the synthetic data. For each inversion, the correlation coeficient between the expected and inverted flow is 
given, the RMS of the expected flow, the expected RMS of the random-noise component, and its measured value. 


may be manipulated by the strength of the noise-regularisation 
term. Generally speaking, the stronger the regularisation, the 
fewer reversals. It is interesting to point out that in the past 
the depth of the supergranulation was estimated from the depth 
where the inv e rted flow reversed its sign IDuvall et al.|([T9 97 ); 
Duvall ( 1998 ); Zhao & Kosovichev ( 20031 ); Svanda etldj i 2009 ) 
reported the reversals depth within a few Mm below the surface. 
Recently, DeGra ve et al. ( |2014b| ) used similar time-distance in¬ 
version as this study to validate it by using the state-of-the-art 
simulation of solar convection and saw the supergranular flow 
reversal, even when it was not present in the simulation. The 
authors claimed that it was due to the shortcomings of the helio- 
seismic inversion methods. 


3.1. Verification of the inversion using synthetic data 

A natural explanation of the issues described above would be 
an error (a programmer’s bug) in the inversion code, or a math¬ 
ematical problem in the inversion itself (e.g., a degeneracy of 
the matrix to be inverted). To eliminate such possibility, a test 
involving the synthetic data was applied, similarly to the tests 


Svanda et al. 


( 2011 ). 


performed by 

A snapshot from the simulation of the convective Sun-like 
flows jUstyugov||2006| was convolved with the appropriate set 
of the sensitivity kernels, in accordance with Eq. ([!]), to create 
a set of synthetic travel-time maps. A random realisation of the 
travel-time noise having the covariance matrix that of the inver¬ 
sion was added to these maps to mimic the random excitation 
of solar waves. In total, three different sets of travel-time maps 
were computed, one for the ridge filters, one for the phase-speed 
filters and one for all filters combined. These travel-time maps 
were inverted by using exactly the same inversion weights as 
used above. 

The resulting flow maps were compared on the pixel-to-pixel 
basis to the flow maps obtained by the smoothing of the simu¬ 
lation with the inversion averaging kernels. Essentially, the left- 
hand side of Eq. was directly compared to the first term of 
the right-hand side of the same equation, denoted as v akern hence¬ 
forth. The root-mean-squared (RMS) value of the random noise 
term (the second term on the right-hand side of Eq. ([2])) was then 
compared to the expected level of the random noise, another out¬ 
put of the inversion. The results are summarised in Table. [I] 

One can see that the inversion behaves as expected. The in¬ 
verted flow maps v inv are highly correlated to the expected flow 
maps v akern . The determined level of the random noise is consis¬ 



Fig. 3. A direct comparison of the expected horizontal flow values to the 
inverted ones on the pixel-to-pixel basis (only each 10th pixel is plotted 
to make the plot simpler). The solid line represents the linear fit to the 
crosses, the thin dotted lines the intervals of the predicted random-noise 
(lcr interval), and the thick dashed line is the line with the unity slope. 



M r+P 

M r 

M p 

Xr+p [S] 

0.24 

0.24 

0.32 

[S] 

0.26 

0.30 

0.33 

xl [s] 

0.20 

0.16 

0.32 


Table 2. Match of the forward-modelled travel times for three consid¬ 
ered models to the observed ones. 


tent with the predicted one from the inversion. The worsts case 
(i.e. depth 6.2 Mm for the phase-speed kernels, which has the 
lowest correlation coefficient and the signal-to-noise ratio) is dis¬ 
played explicitly in Fig. [3] Even this worst case scenario demon¬ 
strates that the inversion behaves well and no sign of a obvious 
bug is found. 

3.2. Does the inverted flow fit the observed travel times? 

In SOLA inversions it is not guaranteed that the inverted flow 
model provides a reasonable fit to the observed travel times. 
It is obvious already from Eq. Q: in order to obtain forward- 
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Fig. 4. Estimated continuous horizontal velocities obtained from four tomographic maps at depths of 1.9, 2.9, 4.3, and 6.4 Mm. The magnitudes 
are saturated at levels +0.5 kms -1 . 


modelled travel times, the knowledge of the continuous flow 
model is needed, which is not returned by SOLA. A proper 
reconstruction is probably possible (a work in progress), how¬ 
ever for a simple illustration a much rougher estimate v® st can be 
made: 

v« st (/, Z) = L l '"' v + ^o)*A"i D (z; Zo), (3) 

Zo 

where = f d 2 r9C^(r,z;zo) is a horizontally averaged av¬ 
eraging kernel. Such an estimate approximately fulfills Eq. 0 
with the random-noise term excluded. The v Q x st components of 
three different models from inversions based on different fil¬ 
tering schemes are shown in Fig. [4] - M P+R for the combined 
ridge+phase-speed filtering scheme, M R for ridge and M P for 
phase-speed filters based inversions respectively. The suspicious 
alterations of the sign of the horizontal flow with depth are 
prominently visible. Models M P+R and M P are similar at depths 
less than 5 Mm and differ by sign deeper down. Model M R is 
significantly different from the other two. 

From these models, forward-modelled travel-time maps 
were computed following Eq. 0 with the noise term neglected. 
Such neglection was fully justified by the ensemble averaging 
technique, which increased the signal-to-noise of the measured 
travel-time maps by more than two orders. A match to the ob¬ 
served travel times r a to the forward-modelled ones rj? m was es¬ 
timated from 

^ = ^Z RMS ( T ? m (r) - Tfl(r) )’ (4) 

a 

where M was the total number of travel-time measurements in¬ 
dexed by a. The larger the worse the forward-modelled 
travel times fit the observed ones. 

The fit was evaluated for each considered vector-flow model 
for all three sets of travel-time measurements separately. The 
results are summarised in Table [2] One would naively expect 
the main-diagonal terms be much smaller than the off-diagonal 
terms. From a visual inspection both the observed and forward- 
modelled travel-time maps it becomes clear that all the fits are 
far from being satisfactory. Dispite being very different (even 
by a sign), all three models fit the observed travel times equally 
(badly). 

The simple reconstruction of the continuous flow field prob¬ 
ably also affects the results, but given the fact that the alteration 
of the sign is visible already in the tomographic maps, the aver¬ 
aging kernels are very confined around the target depth and have 
negligible sidelobes, it can be expected that a more realistic es¬ 
timate of the continuous flow field will not dramatically change 
the conclusions. 
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4. Possible causes 


There are several suspicious contributors possibly responsible 
for the issues. 

Incompatibility of the travel time maps with the inver¬ 
sion weights. As it was pointed out many times in the past by 
various authors, it is crucial that any step in the data process¬ 
ing (mapping, filtering, the way the travel times are measured) 
must be taken into account when computing the corresponding 
sensitivity kernels. In our case, this need was enforced as much 
as possible. The pixel size was exactly the same in case of both 
the data processing and the computation sensitivity kernels, the 
spatio-temporal filtering was done not only using the same code, 
but moreover using the exactly same files, in which the filters 
were stored. The Born sensitivity kernels ([Birch & Gizon|2007| 
are consistent with the |Gizon & Birch| ( |2Q04| ) definition of the 
travel times. The issue with unknown impact are the possible 
non-linearities in both the forward and inverse problems. 

Difference of the power spectrum of the real data and 
the model. The sensitivity kernels are computed using the 
power spectrum, which comes from the model, which is con¬ 
sidered Sun-like ( [Christensen-Dalsgaard et al.|1996[ ). However, 
the eigenfrequencies deviate from the data eigenfrequencies for 
higher order modes (already a mode frequencies are consid¬ 
erably off), the eigenfrequencies are not available for very small 
wave numbers (k < 0.1 Mm -1 ). To investigate the influence of 
this issue, I recomputed sensitivity kernels, travel times, and in¬ 
versions with an additional spatio-temporal filter, which limited 
the power spectrum to the region where the model power spec¬ 
trum matched well to the power spectrum of the used datacubes. 
Hence low wave numbers were filtered out, as was all the sig¬ 
nal of oscillations beyond the ridge. The introduction of this 
additional filter had only minor impact on the results. 

Inaccuracy of sensitivity kernels. Should the sensitivity 
kernels forward-modelled from a reference solar model be dif¬ 
ferent from “those of the real Sun”, we must expect an impact 
of this deviation to the results. The difference is expected to 
be small if the reference model is Sun-like. Unfortunately, there 
does not seem to be a direct way to measure the sensitivity ker¬ 
nels from solar data, perhaps except for the iterative inversions 
( [Hanasoge & Tromp|2014[ H anaso g e|2014| ), introduced to helio¬ 
seismology only recently. There might be the way to at least ver¬ 
ify the total integral of the kernel, which is the work in progress 
and will be reported on in a later paper. This point embodies 
both the effects of the differences between the reference model 
and the real Sun and the possible issues with the kernel computa¬ 
tion, both theoretical and numerical. This point is also supported 
by peGrave et al.| ( |2014b| ), where the authors using a realistic 
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Fig. 5. The contributions of various filtering schemes to the inversion averaging kernel in case of the combined filtering scheme. 


convection simulation noticed a significantly decreasing corre¬ 
lation between the measured and forward-modelled travel times 
for ridge filters with higher orders ( p 3 and beyond) and also for 
larger phase speeds. 

Mathematically posed problem. This possibility is closely 
related to the previous point. Let’s say that the sensitivity kernels 
are accurate (again, by “accurate” it is ment “those of the real 
Sun”) to within 1%. The combined inversion combines 405 such 
kernels. These 1% errors translate through the inversion and may 
be amplified to an unknown extent. There are some hints that this 
is likely the most probably cause. The contributions of the aver¬ 
aging kernel from two different filtering schemes are displayed 
in Fig. [5] One can see that in all four cases discussed in this re¬ 
search note the two components from different filtering schemes 
largely subtract from each other. Should one filtering scheme 
have an unknown systematical bias, the resulting “real” aver¬ 
aging kernel (i.e. the averaging kernels obtained by convolving 
inversion weights with the sensitivity kernels “of the real Sun”, 
which are unknown) is different and that would easily explain 
even the change of sign of the flow. A curious reader may object 
that the way out is not to combine the two filtering schemes. It 
is shown in Fig. [I] that the inverted horizontal flow based solely 
on the phase-speed filtering scheme reverses its sign between the 


depths of 1.9 Mm and 2.9 Mm, which is not realistic. One has to 
bear in mind that the subtractions occur in every inversion, e.g. 
the contribution of point-to-annulus measurements is subtracted 
from the contribution of the point-to-quadrant measurements etc. 
Another indication for this explanation is that the magnitude of 
the inverted flow which uses the phase-speed filters seems to be 
a bit unrealistic (around 2 km s -1 for the depth of 1.9 Mm). 

Both last points implicate that the “real” inversion averaging 
kernels may be very different from those predicted by the in¬ 
version, perhaps even having a negative total integral. Then also 
the cross-talk contribution is not constrained, which is estimated 
from the cross-talk components of the averaging kernel. 


5. Lessons learned 


This work is based on many tens of thousands of CPU-hours 
trial-and-error runs. In case of the inver sions suitable for the to¬ 
mography of the flow snapshot (see e.g. 


Svanda 


2013) the issues 


are not prominently visible, however the question is, whether the 
stronger regularisation of the solution (about the random-noise 
term in this case) removes the issues or rather hides them. By 
studying the literature on flow inversions in supergranules one 
has to conclude that all the inversions except for the inversion 
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of |Woodard] ( [20071 ) indicated the reversal of the horizontal flow 
at various depths. It may be seen suspicious, because the state- 
of-the-art simulations ( [Stein et al.|2009[ |Rempel et al.|2009| ) do 
not indicate such reversals. It may also easily be that the large- 
amplitude flows in supergranulati on recently reported by time- 
distance inversions (Svanda 2012) are artefacts of similar issues, 
as in that case / mode and p\ and P 2 ridges were used, necessar¬ 
ily leading to subtractions in the inversion. 

From comparisons of the inversion results with the direct sur¬ 
face measurements (e.g. | Ambroz|2005l|Georgobiani et al.|2007 


Svanda et al. 


2007, Svanda et al.||2013) it seems that the very 


near-surface inversions (hence involving / modes or acoustic 
waves with small phase speeds) do not suffer from the discussed 
problems. As we showed recently (Svanda et al. 2013), the flow 
inversion using the / mode ridge is not only highly correlated 
with the inferences from the surface granule tracking, but it also 
provides the properly scaled magnitude of the flow. Such vali¬ 
dations against the independently obtained measurements fully 
justify the scientific results obtained on surface flow fields. It 
has to be pointed out that no larger subtractions occured in those 
inversions. 


Acknowledgements. This work was supported by the Czech Science Founda¬ 
tion (grant P209/12/P568). All computations were performed using the Sun- 
quake compute cluster at Astronomical Institute of Academy of Sciences in 
Ondrejov, the tracked and mapped datacubes were obtained from data processing 
pipelines at Max-Planck-Institut flir Sonnensystemforschung (MPS), Gottingen, 
Germany, which is funded by the German Aerospace Center (DLR). The solar 
measurements were kindly provided by the HMI consortium. The point-to-point 
travel-time sensitivity kernels were obtained using the code of Aaron Birch de¬ 
ployed at the MPS. This work was done within the institute research project 
RVO:67985815 to Astronomical Institute of Czech Academy of Sciences. This 
work benefits from couloir discussions following the talk Svanda: Current issues 
with time-distance inversions for flows at As U at LWS Helioseismology Work¬ 
shop #4: Solar Subsurface Flows from Helioseismology: Problems and Prospects 
which was held at Stanford University, July 21-23, 2014. 


Stein, R. F., Nordlund, A., Georgoviani, D., Benson, D., & Schaffenberger, W. 
2009, in Astronomical Society of the Pacific Conference Series, Vol. 416, 
Solar-Stellar Dynamos as Revealed by Helio- and Asteroseismology: GONG 
2008/S OHO 21, ed. M. Dikpati, T. Arentoft, I. Gonzalez Hernandez, C. Lind¬ 
sey, & F. Hill, 421 

Svanda, M. 2012, ApJ, 759, L29 

Svanda, M., Gizon, L., Hanasoge, S. M., & Ustyugov, S. D. 2011, A&A, 530, 
A148 

Svanda, M., Klvana, M., Sobotka, M., Kosovichev, A. G., & Duvall, T. L. 2009, 
New Astron., 14, 429 

Svanda, M., Zhao, J„ & Kosovichev, A. G. 2007, Sol. Phys., 241, 27 

Ustyugov, S. D. 2006, in Astronomical Society of the Pacific Conference Series, 
Vol. 354, Solar MHD Theory and Observations: A High Spatial Resolution 
Perspective, ed. J. Leibacher, R. F. Stein, & H. Uitenbroek, 115 

Svanda, M. 2013, ApJ, 775, 7 

Svanda, M., Roudier, T., Rieutord, M., Burston, R., & Gizon, L. 2013, ApJ, 771, 
32 

Svanda, M., Sobotka, M., & Barta, T. 2014, ApJ, 790, 135 

Woodard, M. F. 2007, ApJ, 668, 1189 

Zhao, J. & Kosovichev, A. G. 2003, in ESA Special Publication, Vol. 517, 
GONG+ 2002. Local and Global Helioseismology: the Present and Future, 
ed. H. Sawaya-Lacoste, 417-420 

Zhao, J., Kosovichev, A. G., & Duvall, Jr., T. L. 2001, ApJ, 557, 384 


References 

Ambroz, P. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 
346, Large-scale Structures and their Role in Solar Activity, ed. K. Sankara- 
subramanian, M. Penn, & A. Pevtsov, 3 

Birch, A. C., Braun, D. C., Leka, K. D., Barnes, G., & Javornik, B. 2013, ApJ, 
762, 131 

Birch, A. C. & Gizon, L. 2007, Astron. Nachr., 328, 228 
Christensen-Dalsgaard, J., Dappen, W., Ajukov, S. V., et al. 1996, Science, 272, 
1286 

Couvidat, S., Birch, A. C., & Kosovichev, A. G. 2006, ApJ, 640, 516 
Couvidat, S., Gizon, L., Birch, A. C., Larsen, R. M., & Kosovichev, A. G. 2005, 
ApJS, 158, 217 

DeGrave, K., Jackiewicz, J., & Rempel, M. 2014a, ApJ, 794, 18 
DeGrave, K., Jackiewicz, J., & Rempel, M. 2014b, ApJ, 788, 127 
Duvall, Jr., T. L. 1998, in ESA Special Publication, Vol. 418, Structure and Dy¬ 
namics of the Interior of the Sun and Sun-like Stars, ed. S. Korzennik, 581 
Duvall, Jr., T. L. & Birch, A. C. 2010, ApJ, 725, L47 
Duvall, Jr., T. L., Birch, A. C., & Gizon, L. 2006, ApJ, 646, 553 
Duvall, Jr., T. L., Jefferies, S. M., Harvey, J. W., & Pomerantz, M. A. 1993, 
Nature, 362, 430 

Duvall, Jr., T. L., Kosovichev, A. G., Scherrer, P. H., et al. 1997, Sol. Phys., 170, 
63 

Georgobiani, D., Zhao, J., Kosovichev, A. G., et al. 2007, ApJ, 657, 1157 
Gizon, L. & Birch, A. C. 2004, ApJ, 614, 472 
Gizon, L., Birch, A. C., & Spruit, H. C. 2010, ARA&A, 48, 289 
Hanasoge, S. M. 2014, ArXiv e-prints 
Hanasoge, S. M. & Tromp, J. 2014, ApJ, 784, 69 
Jackiewicz, J., Birch, A. C., Gizon, L., et al. 2012, Sol. Phys., 276, 19 
Jackiewicz, J., Gizon, L., & Birch, A. C. 2008, Sol. Phys., 251, 381 
Kosovichev, A. G. & Duvall, Jr., T. L. 1997, in Astrophysics and Space Science 
Library, Vol. 225, SCORe’96 : Solar Convection and Oscillations and their 
Relationship, ed. F. P. Pijpers, J. Christensen-Dalsgaard, & C. S. Rosenthal, 
241-260 

Rempel, M., Schiissler, M., & Knolker, M. 2009, ApJ, 691, 640 
Article number, page 8 of[8] 






















